Spin self-rephasing and very long coherence times in a trapped atomic ensemble 
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We perform Ramsey spectroscopy on the ground state of ultra-cold ^^Rb atoms magnetically trapped on 
a chip in the Knudsen regime. Field inhomogeneities over the sample should limit the 1 /e contrast decay 
time to about 3 s, while decay times of 58 ± 12 s are actually observed. We explain this surprising result 
by a spin self-rephasing mechanism induced by the identical spin rotation effect originating from particle 
indistinguishability. We propose a theory of this synchronization mechanism and obtain good agreement 
with the experimental observations. The effect is general and may appear in other physical systems. 



In atomic clocks and otherprecision techniques based on 
atomic spin manipulation 111], a central requirement is to 
preserve the coherence of a state superposition over long 
times. Understanding how coherence decays in a given 
system is important for these applications, and is a touch- 
stone of understanding its dynamics. In trapped ensem- 
bles, an inhomogeneous shift A(r) of the transition fre- 
quency occurs due to the trapping potential and to atomic 
interactions. Different atoms explore different regions of 
this shift landscape, and so their spins precess at different 
rates. This leads to dephasing at a rate determined by the 
characteristic inhomogeneity Aq of A(r) over the ensem- 
ble. Various mechanisms have been exploited to reduce this 
dephasing. Examples are "magic fields" that strongly re- 
duce the field dependence for a specific transition |y,[3|], or 
the mutual compensation scheme successfully employed in 
ultracold ®^Rb [41], where the trap-induced inhomogeneity 
can be adjusted to nearly cancel the collisional mean-field 
inhomogeneity. All such mechanisms however, including 
the motional narrowing well known in nuclear magnetic 
resonance, have in common that the dephasing is merely 
slowed down, but never reversed, and the transverse polar- 
ization remains a steadily decreasing function of time. 

Here we present measurements on a trapped ensemble 
of ^^Rb atoms with two internal levels equivalent to a spin 
1/2. Atomic interactions cause a spontaneous re-phasing 
of the spins, observed as a much longer decay time and re- 
vivals of Ramsey contrast. We are also able to extend the 
coherence time by more than an order of magnitude beyond 
the 2 to 3 s previously achieved on this system Oilsl]. We 
explain these remarkable results by a very general mecha- 
nism based on the identical spin rotation effect (ISRE) that 
occurs during collisions in the forward direction between 
two identical particles |@] - an equivalent description can 
be given in terms of the exchange mean-field experienced 
by the atoms ||7|]. This effect is known to cause transient 



(c) 



spin waves 13,11-14], a deleterious phenomenon if one is 
interested in long coherence times. In contrast to those ex- 
periments however, we are working in a regime where both 
(i) the ISRE rate (exchange rate) lo^x/'^t^ = 2?i|aoi \n/m is 




FIG. 1: Two classes of atoms (red and blue) precess at differ- 
ent rates. Their Bloch vectors were initially parallel, but have 
started to dephase (a). The ISRE then makes both vectors rotate 
around their sum (grey) (b). When this rotation reaches tt, the 
fast-precessing spin (red) lags behind the slow one (blue), which 
tends to rephase them as in a spin echo (c). 



larger than the inhomogeneity Aq and in) the rate of lateral 
elastic collisions 7^ = {"ilyjlr /'i)a^^nvx is much lower 
than both the trap frequencies (Knudsen regime) and ujg^. 
Here, Oqi is the relevant scattering length ||24|| . n the aver- 
age density and Vt = yJkoT jm the thermal velocity for 
atomic mass m. In this regime, we find that the ISRE intro- 
duces efficient synchronization between atoms with differ- 
ent spin precession rates. Indeed, solving a kinetic equation 
for the spin variables based on the ISRE, we obtain good 
agreement with the data. These findings are reminiscent of 
earlier calculations for a trapped gas which predict local- 
ized polarization revivals [12] and synchronization within 
spatial domains [ICj] in a different (hydrodynamic) regime. 



We start with a simple model (figure [T]), where as in lll3h 
we divide the atoms into two classes having fast and slow 
transverse spin precession rates. Whether an atom is in 
the fast or slow class depends on the average of A(r) it 
experiences. In the second experiment below, for example, 
this depends on its orbital energy. If the class-changing 
events are rare (7^ < a;^,^), an atom remains in its class for 
a long time and the two classes start to dephase. The effect 
of the ISRE is then simply to make the spin polarizations 
of the slow and fast class turn around their sum, reversing 
their phase differences on the timescale of Tr/wgr. When 
this rotation reaches vr, the two transverse polarizations are 
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exchanged, while each of them continues to precess at the 
same rate as before, so that they start to rephase. 

The ISRE does not change the sum of the polarizations: 
it just reverts the correlation between the spin directions 
and their precession rate. This initiates the rephasing pro- 
cess that increases the transverse polarization. Of course, 
in reality dephasing and refocussing occur simultaneously; 
the result is a synchronization mechanism that is analogous 
to a negative feedback effect. It is also reminiscent of a 
spin echo, but here the rephasing of the spins is due to an 
internal effect instead of externally applied RF pulses. 

For a self-rephasing regime to occur, the ISRE must 
revert the spins before lateral collisions cause atoms to 
change class - that is, w^v/vr > 7c - and before the de- 
phasing reaches vr, oj^v > Aq (otherwise it accelerates the 
dephasing). If S> Aq, any dephasing will be immedi- 
ately refocussed, leading to tight synchronization and long 
coherence times; if > Aq, a significant phase spread 
can occur before the spins rephase, leading to loss and re- 
vival of the total polarization. Performing two experiments, 
we have indeed observed these behaviours, and found good 
agreement with our quantitative calculations. 

We perform Ramsey spectroscopy on the |0) = 
l-F = l,mp = —1) to |1) = |F = 2, nip = 1) hyperfine 
transition of * ''Rb atoms magnetically trapped on a chip. 
The differential Zeeman shift Az{t) is quadratic in mag- 
netic field around a "magic" value Bm = 3.228917(3) G 
iH] and quartic in position. The s-wave scattering lengths 
Ooo, flu and agi differ by less than 5% lITsIl , so that the 
mean-field shift of the transition frequency is already small 
with A,„/(r)/27r = -0.4 Hz x n{r)/10^^cm-^ for an 
equal superposition of |0) and where n(r) is the local 
density. Additionally, we employ the mutual compensa- 
tion technique Jsl]: adjusting the magnetic field at the trap 
center Bq to a value slightly below Bm leads to a standard 
deviation of A(r) = Az{r) + Aj„j(r) averaged over the 
atomic cloud that is of order Aq = 27r x 0.08 Hz lfl6ll . 

The experiment is designed with the goal of a "trapped- 
atom clock on a chip" (TACC), aiming at a stability in the 
lower 10~^'^s~^/^ range. The setup is described in detail 
in ifT/ll . and is similar to that of ifisll . It incorporates a 
two-layer atom chip with a coplanar waveguide (CPW). 
Atom preparation involves the usual laser and evapora- 
tive cooling steps ifigll . Within less than 10s, we obtain 
A'at = 5 X 10^ to 10^ trapped atoms in state |0) at a 
temperature of T = 175(6) nK, which is at least 30 nK 
above the onset of Bose Einstein condensation (BEC) for 
our densities. We use a two-photon, microwave (MW) and 
radiofrequency (RF) excitation to drive the clock transition. 
The MW frequency is detuned ~ 500 kHz above the 1 0) to 
\F = 2, rup = 0) transition. The detuning from the two- 
photon resonance is set to Aj^/27r = 3.6 Hz. The MW 
is generated by a home-built synthesizer having very low 
phase noise ||20|| ; a commercial direct digital synthesizer 
provides the RF. Both are locked to a hydrogen maser of 
relative frequency stability lO^^^r"^/^ up to 1000s ll2lll . 



Both signals are injected into structures on the chip with 
powers ~ dBm (MW) ~5 dBm (RF). Trapped below the 
CPW at a distance Zq, the atoms interact with its evanes- 
cent field and perform two-photon Rabi oscillations at a 
rate 0/^. In each cycle, we detect the number of atoms 
in both clock states, A'^i and Nq, by absorption imaging 
without and with repump light after a time of flight of 9 ms 
and 13 ms respectively. We use absorption imaging close 
to saturation intensity ||23]. Intensity and magnification are 
carefully calibrated, and we have verified that cross talk be- 
tween the detected and A'^i is negligible. We estimate 
the calibration error on the absolute number of atoms to 
5%. The statistical error is 140 atoms per state. 

In a first experiment, the magnetic trap has frequencies 
{uj^,Uy,u^]/2'K = {32(1), 97.5(2.5), 121(1)} Hz; for 
this trap = 156 /im and Vtji = 2tt x 164 Hz. Bq is op- 
timized roughly by maximizing the fringe contrast after a 
Ramsey time Tr = 2 s. This leads to Bq = 3.1626(7) G, 
as measured by RF-induced atom loss on a small BEC 
(^at ~ 3000). The two 7r/2 pulses, spaced by a vari- 
able Ramsey time Tr, are applied at the end of a period 
of fixed length of Tp = 5.02 s, and preceded by a hold- 
ing time Th = Tp — during which the atoms are in 
|0). The transition probability Ni/{Nq + iVi) is plotted 
in figure |2a) as a function of Tr. The contrast at Tp = 
is 89%, probably limited by a small inhomogeneity of the 
Rabi frequency. At Tp = 5 s the contrast still remains 
above 82%. Assuming exponential decay, the 1/e time is 
58 lb 12 s, much longer than the 2.75 s predicted from the 
mutual compensation scheme liT6ll . The true dephasing is 
even slower, since the atom loss rates 70 and 71 of |0) and 
|1) differ (fig. |2b)); asymmetric loss leads to a decay of 
Ramsey contrast that is independent of dephasing. 

The very long coherence time can be understood as a 
tight synchronization due to self-rephasing. The density 
n ~ 1 (from now on in units of lO^^cm"'^) gives an ISRE 
rate ti;gv/2vr ~ 8 Hz, much larger than the inhomogene- 
ity Ao/27r ~ 0.08 Hz, so that the fast and slow spins are 
swapped when the dephasing is still small. Moreover, since 
7c ~ 2s^^, the correlations between precession rate and 
accumulated dephasing remain intact over long times. 

In a second experiment, we deliberately increase 
Aq by detuning Bq away from the optimal value. 
We choose Bq = 3.7562(9) G, such that the in- 
homogeneity of the transverse precession rate is now 
well approximated by a parabolic spatial dependence 
A(r) « Aoiix/xpf + {y/yTy + {z/zTy) with 
{xp,yp,ZT} = VT/{uj^,u}y,UJ,}. Here, Ao = 
(x|/2)(5^A) « 27r (1.2 + 0.1 X n) Hzby averaging over 
the gaussian density profile n{r). Note that Aq is also lin- 
ear in temperature, which is constant in our experiment. 

The new Bq leads to trap frequencies 
{uJa^,L0y,L0,}/27r = {31.30(5), 92.0(5), 117.0(3)} Hz; 
Zq = 151 //m and Qp, = 2tt x 83 Hz. The total time 
is reduced to = Is and the Ramsey time is varied 
between 15 and 500 ms. For each Tp, we vary the detuning 
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FIG. 2: First experiment: (a) Ramsey fringes in the time domain for Aq jl-n ^0.08 Hz. The normalized transition probablity is plotted 
as a function of the Ramsey free evolution time. The initial contrast is 89%, probably limited by an inhomogeneous Rabi frequency. At 
Tfj = 5 s, the contrast is 82% resulting in a 1/e time of 58± 12 s. (b) Total number of detected atoms. Although the total time is kept 
constant, this number depends on Tfj because the loss rates of |0) and |1) differ. The data is fitted by A'o + A^i = N^/'i (1 + e~"^"/^) 
with Nt = 24.8(2) x 10^ atoms and r = 8.7(2) s. 



in 30 steps and extract the fringe contrast. We repeat 
the measurement for different n spanning one order of 
magnitude. To vary n, we vary the atom number by 
changing the MOT loading time. (We have checked that 
the temperature is independent of n within our measure- 
ment precision.) The lowest density is n = 0.2, which 
corresponds to u}^^/2t: = 1.5 Hz and 7^, = 0.4s~^. 

The results are plotted in figure [3] All datasets show 
initial rapid contrast decay between and 50 ms. For n = 
0.2, it vanishes completely at ~ 200 ms. Increasing the 
density stabilizes the contrast beyond 200 ms and leads to 
revivals. The contrast improves with increasing density. 
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FIG. 3: Second experiment: Ramsey fringe contrast as 
a function of Tfj for increased inhomogeneity Ao/27r f» 
1.2 Hz. From bottom to top, colors correspond to densi- 
ties n = {0.2,0.8,1.1,1.9,2.6} x lO^^ cm"'"' (or Nat = 
{7.9, 28, 41, 70, 95} x 10^). For each n, the data is normalized 
to the initial contrast at = 0. Solid lines are numerical so- 
lutions of the kinetic equation for Ao/27r ?s 2 Hz taking into 
account the ISRE {ujex/'i''^ ~ 4.5Hz x n) and lateral collisions 
(7c ~ 2.1s^^ X n). The red dotted line is a calculation for n — 
2.6 (same density as the red data points) with lateral collisions 
(7c = 5s~^), but without ISRE. The red dashed line is a calcula- 
tion without lateral collisions but with ISRE (tje,i/27r — 12 Hz). 



indicating an interaction effect. The time of the first revival 
depends on density with Trevivai ~ —0.02 s -|- 0.3 s/n. This 
agrees to within a factor of 2 with the estimation 211 juje^ = 
0.13s/n from our simple model above. 

To interpret this data, we now perform a quantitative cal- 
culation. As in ifTll [l2ll . the motion in the trap is treated 
semiclassically, while a full quantum treatment of the spin 
variables is included in a kinetic equation for a density op- 
erator y6(r, p, t), characterized by the usual Bloch vector 
S = Tr{pa/2). The vector components of S and the 
Pauli matrices a point along the three unit vectors of the 
Bloch sphere {u^i, Uj^2i U||}. To describe the harmonic 
oscillator in phase space, instead of the usual {x,Px} 
(and similarly for y and z), we use energy-angle variables 
{E^ = (x^ + pl)/2,a,j: = arctan(pj./x)}, where E^, 
X and are respectively in units of ksT, Xt and tuvt- 
As the oscillatory motion is fast compared to the spin dy- 
namics (oJx.y.z ^ '^e.vi 7c)> we average S(E,Q!,t) 
over the angles cx = {a:c,Oiy,az) and obtain an equa- 
tion for a spin density S(E, t) in energy space only, where 
E = (E^, Ey, E^) Jq^ . This description is inspired by 
that developed in [g] for a quasi- ID gas in the Knudsen 
regime, but generalized to 3D and with the effect of lat- 
eral collisions giving rise to a damping term oc 7c. At 
this stage, the equation still depends on each Ei separately 
through the frequency shift A(E). However, since our 
A(r) is parabolic in r, it becomes A(E) = AqE when 
averaged over the angles and only depends on the total en- 
ergy E = Ej. + Ey + Ez- This allows us to write an 
equation for an energy-isotropic spin density S{E, t): 



dtS{E,t) + j4S{E,t)-S{t)] ^ 
dE'—e-'''K{E,E')S{E',t) 



A(i?)u|| 
X S(^,t)(l) 



where S = dE^e ^S{E) is the average spin and 
E'^/2 is the 3D harmonic oscillator density of states. 
Upon angle-averaging, the spin mean-field, which is lo- 
cal in position, becomes long-ranged in energy space with 
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a kernel K{E,E'). In ID, this kernel Ki{E,E') ^ 
[max{E,E')\E - ||^. \n 3D, it is also long- 

ranged but cumbersome to use. For simplicity, when solv- 
ing ([ll numerically, we approximate ll26ll it as infinite- 
ranged K(E, E') ~ 1. Eq. ([B does not include atom loss. 
It contains three experimentally tunable parameters: the in- 
homogeneity Ao(T, n, ...), the ISRE rate uje^ oc |aoi|n 
and the lateral collision rate 7c oc a^^n^/T. 

To compare this theory to the experiment, we first con- 
sider the lowest-density dataset (n = 0.2), where both the 
ISRE and lateral collisions are negligible (Aq ^ oOex, 7c)- 
When they vanish, the kinetic equation is easily solved an- 
alytically. For an initial condition S(i?, 0) = uj^i, the 
contrast is |S(t)| = (1 + (Aot)^)"^/^. Fitting the exper- 
imental contrast with this result, we find, Ao/27r ~ 2 Hz, 
not far from its expected value Ao/27r w 1.2 Hz. We 
then numerically solve eq. ([Tl for the other densities 
n = {0.8, 1.1, 1.9,2.6}. Choosing Ao/27r within 5% of 
the above 2 Hz, 7c = {?,2^/^T /2>)a1'^nvT ~ 2.1s^^ x n as 
predicted and w„y27r ~ 0.6 x 2h\a^i\n/m ~ 4.5Hz x n 
reproduces all data well (solid lines in fig. |3]l. The renor- 
malization of the ISRE rate by a factor 0.6 results from 
the overestimation of the synchronization effect through 
the infinite -range approximation. 

If we set ujex = 0, the theory predicts short coherence 
time and no revivals (dotted line in fig.O for all densities, 
confirming that the ISRE is responsible for the revivals. 
We also solve the n = 2.6 case without lateral collisions, 
7c = 0: the dashed line in fig. |3] shows an initial drop and 
revival, but continues to oscillate around a constant value. 
The lateral collisions are therefore responsible for the slow 
decay of the contrast at long times. We note that, in contrast 
to the mechanism of ifioll . here the spin synchronization 
does not result from a simple compensation of the inhomo- 
geneous longitudinal field A(£')u|| by the exchange mean 
field uje,x^{t) as they are orthogonal. Also, while Ref. ll23ll 
discusses the effects of the inhomogeneities of the probe 
field during the 7r/2 pulses, our work deals with the effects 
of the static field between them (free precession). 

In conclusion, the observed spin self-rephasing and syn- 
chronization are relatively robust when the inhomogeneity 
to compensate is not too large: in the first experiment, a 
large orientation is still present after more than 10 velocity- 
changing collisions. The effect may occur in any sample 
where the exchange rate is larger than the characteris- 
tic dephasing rate Aq and the lateral collision rate 7c. The 
ratio between these rates is tunable through the temperature 
dependence of w^v/ 7c and through Feshbach resonances. It 
should be interesting to investigate whether the mechanism 
occurs in systems such as optical lattice clocks jH. 
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